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Abstract 

Due to their highly random nature, vegetation canopies can be modeled using the incoherent transport theory 
for active and passive remote sensing applications. Agricultural vegetation canopies however are generally more 
structured than natural vegetation. The inherent row structure in agricultural canopies induces coherence effects 
disregarded by the transport theory. The objective of this study is to demonstrate, via Monte-Carlo simulations, these 
coherence effects on L-band scattering and thermal emission from corn canopies consisting of only stalks. 

1. Introduction 


In this study, we demonstrate via Monte-Carlo (MC) simulations, the effect coherence induced by vegetation 
canopy row structure on active and passive remote sensing of the canopy. For this purpose, a simple defoliated corn 
canopy consisting of only stalks is considered. In order to imitate the canopy row structure, the stalks are distributed 
with a quasi-periodic distribution (QPD) where they are first placed on a periodic lattice and later given random 
displacements from their positions on the periodic lattice. Monte-Carlo simulations are performed for the mean 
backscattering coefficient and brightness temperature at F-band for sensors looking along the row direction and 
perpendicular to the row direction to observe the effect of coherence. Results are also compared with those from 
uniform canopies with stalks distributed uniformly throughout the canopy. Discussion on the implications of the 
effect of coherence is given at the end. 

2. Quasi-Periodic Canopy Model 

The quasi-periodic canopy model consists of stalks (modeled as finite length dielectric cylinders) randomly 
placed with a quasi-periodic distribution over a perfect electric conductor (PEC) plane. For simplicity, the leaves 
and cobs are assumed to be removed from the stalks. The QPD can be regarded as a non-ideal periodic structure. It 
is obtained from the periodic lattice of the cylinders by randomly displacing the cylinders from their ideal positions 
on the periodic lattice. The cylinders are displaced within a bounding region according to a prescribed distribution 
with zero mean. In this manner, all cylinders are identically distributed about their ideal positions on the lattice. Note 
that the periodic canopy can be obtained from the quasi-periodic canopy by choosing a Dirac delta distribution. A 
top view of the quasi-periodic canopy with circular bounding regions of radius b , is depicted in Fig.l. In the figure, 
A x and A v correspond to the row spacing and the in-row cylinder spacing respectively. The canopy composed of 
finite length cylinders over a PEC plane is depicted in Fig.2. 


3. Analysis 


For the active problem, we consider a unit amplitude plane wave with wavenumber, /c 0 , incident on the canopy 
in the direction defined by the angle pair (6 it 0j) where the angle Q t is measured from the z-axis perpendicular to the 
PEC plane and the azimuth angle (p t is measured from x-axis in the row direction. The corresponding incidence vector 
k is shown in Fig 2. The mean bistatic cross section of the canopy is defined as: 
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Figure 1. Top view of a quasi-periodic canopy. Figure 2. Model for corn canopy composed of finite length 

dielectric cylinders over a PEC plane. 


In the above, r is the distance of observation from the canopy, E op and E sq are the /^-polarized incident and q- 
polarized scattered field amplitudes respectively (; p , q E { h , v }), A is the area of the canopy and ( ) denotes the 
average over cylinder configurations. The term E S f q is the ^-polarized fluctuating part of the scattered field defined 
as E S f >q = E s q — ( E s q ). The quantity of interest in this study is the mean backscattering coefficient obtained from 
(1) with k! = —k. Considering the mean scattered field, ( E s q ), it is clear that a QPD canopy of infinite extent is 
invariant under translations by multiples of A x or A y in the x and y directions respectively. This will lead to space 
harmonics for the mean scattered field determined by the angles 
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where m,n E Z i.e. integer. For a canopy of finite extent, the space harmonics will persist but with finite angular 
width inversely related to the canopy extent. 

For the passive problem, the emissivity of the canopy is given by the generalized Kirchhoff law i.e. the 
emissivity of the canopy is equal to its absorbtivity under thermal equilibrium. Under thermal equilibrium the active 
and passive problems are related as follows: emission in the ( 6 if (p t ) direction is obtained from the absorption by the 
stalks due to a unit amplitude plane wave incident on the canopy in the (6 if 0*) direction. Precisely, the /^-polarized 
emissivity of the canopy is given by: 
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where (J° ap (Qi, 0;) is the absorbtivity of the canopy due to the unit amplitude /^-polarized plane wave incident in the 
(6 if 0 j) direction and d£l = sin Q dOdcp is the differential solid angle. The above expression is commonly known as 
Peake’s principle [1]. Note that since the cylinders (stalks) are placed over a PEC plane, absorption (hence emission) 
will correspond only to power absorbed (emitted) by the cylinders. 

In order to compute the scattered and absorbed powers, the fields on the surface of the cylinders must be 
known. For this purpose, using a cylindrical harmonic expansion, a set of self-consistent (Foldy-Lax) equations is 
constructed for the coefficients of the surface fields of the cylinders. We emphasize that all multiple scattering is 
taken into account in this approach. The solution to the self-consistent equations requires a large matrix inversion 
that is performed using a bi-conjugate gradient method in conjunction with an FFT algorithm to speed up matrix- 
vector multiplication. Details of the implementation can be found in [2] and [3]. Once the surface fields are 
obtained, the backscattering coefficient and the emissivity can be found using (1) and (3) respectively. In this study, 
emissivity is alternatively obtained from the power absorbed by the cylinders via a surface integral of the Poynting 
vector on the surface of the cylinders [3]. 

The averaging of the backscattering coefficient and emissivity is performed via MC averaging. Each 
realization of the simulation procedure consists of (i) generating a sample quasi-periodic canopy, (ii) solving for the 
coefficients of the cylinder surface fields (iii) computing the backscattering coefficient and the absorption cross 
section of the canopy. Mean backscattering and emissivity are obtained by averaging over all MC realizations. 



4. Simulation Results and Discussion 


In this section, simulation results are given for the mean backscattering coefficients and mean emissions for 
two quasi-periodic com canopies. Simulations are performed for both parallel-to-row and perpendicular-to-row look 
directions. In addition, we associate a uniform canopy with each quasi-periodic canopy such that the canopy has the 
same density as its associated quasi-periodic canopy. The stalks in the canopy have a uniform identical and 
independent distribution over the canopy area, A. Simulations are performed for these uniform canopies to facilitate 
comparison with their associated quasi-periodic canopies. 

The simulations are performed at 1.5 GHz for vertical polarization. Each MC simulation is the result of 
averaging over 50 realizations. The dielectric cylinders representing the stalks have radius a = 1cm, length L = 1 m 
and relative permittivity 8 r = 15 + iS and are identical for all canopies. A total of 512 cylinders are used in each 
canopy. Other relevant canopy parameters are given in Table 1. 


Table 1. Canopy Parameters forMC Simulations 



Canopy- 1 

Canopy-2 

Stalk Density #/m 2 

8 

6.7 

Row spacing (m), A x 

0.50 

0.75 

Average stalk spacing (m), A v 

0.25 

0.20 

Bounding region radius, b (m) 

0.10 

0.08 


Simulations results of mean backscattering coefficient vs. incidence angle for Canopy- 1 and Canopy-2 are 
given in Fig. 3a and Fig. 3b respectively. In both figures three curves are plotted. The curves denoted by QPD(||) and 
QPD(±) correspond to quasi-periodic canopy results with parallel-to-row and perpendicular-to-row looks at the 
canopy respectively. The curve denoted by UIID corresponds to results of the associated uniform canopy. It is seen 
from Fig. 3a that for Canopy- 1 the three curves are within 2 dB of each other. For Canopy-2 however, above 30° 
incidence there is about 5 dB difference between backscattering coefficients for the perpendicular look at the QPD 
canopy and the uniform canopy. The parallel look at the QPD canopy, on the other hand, has a conspicuous peak at 
30° incidence angle. This peak corresponds to a space harmonic occurring in the backscattering direction as can be 
computed from (2) with values from Table 1 corresponding to Canopy-2. 

Simulation results of mean brightness temperature vs. incidence (observation) angle for Canopy- 1 and 
Canopy-2 are given in Fig. 4a and Fig. 4b respectively. The curves in Fig. 4a show that the brightness temperatures 
for the parallel and perpendicular looks at Canopy- 1 are very similar over the complete angle range with a maximum 
difference of about 10 K. The brightness temperatures for both parallel and perpendicular looks at the QPD canopy 
differ from those of the uniform canopy with a maximum difference at 6 t = 40° of about 36 K between the parallel 
look and the uniform canopy and a maximum difference of about 26 K between the perpendicular look and the 




Figure 3. MC simulations of backscattering coefficient vs. incidence angle: (a) results for Canopy- 1 
(b) results for Canopy-2. 




Figure 4. MC simulations of brightness temperature vs. observation angle: (a) results for Canopy- 1 
(b) results for Canopy-2. 

uniform canopy. For Canopy-2, Fig. 4b shows that the brightness temperatures for the parallel and perpendicular 
looks are close up to 0 t = 30°. As the incidence angle is increased, a gradual increasing in difference is observed. 
The maximum difference in the brightness temperatures between these looks is approximately 63 K at 6 t = 50°. It is 
also observed that the brightness temperature of the uniform canopy is close to that of the QPD canopy up to 6 t = 
30° and at higher angles falls between the curves of the parallel and perpendicular looks at the QPD canopy. 

5. Conclusion 

Simulation results confirm the effect of coherence in both active and passive remote sensing of quasi-periodic 
canopies. In the active problem, backscattering can be substantially influenced by coherence if a space harmonic 
coincides with the backscattering direction in which case the uniform canopy assumption fails. This implies that for 
such a configuration, the incoherent transport theory analysis, where a uniform canopy is generally assumed, will 
not be adequate. An analysis using transport theory with non-uniform canopy assumption is given in [4] where it is 
also reported that taking coherence into account yields better results. 

In the passive problem, the effect of coherence is generally visible at higher observation angles. The results 
also render transport theory analysis of emission from such quasi-periodic questionable. It is also worth noting, a 
direct relation between emission and space harmonics is not deduced, nevertheless they enjoy an indirect 
relationship through Peake’s principle given by (3). 

The above results are obtained for defoliated corn canopies. In a real canopy, the leaves will attenuate 
multiple scattering between stalks and hence decrease the effect of coherence and consequently the dependence of 
emission on the row structure. The incorporation of the leaves into the current model requires further study. 
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